geno<-read.table("outfile",sep="\t")
geno<-read.table("sample_outfile",sep="\t")
geno<-geno[,2]
expr<-read.table("sample_outfile_expr",sep="\t")
expr<-expr[,2]
z<- 1
while(z<=length(expr))
{
kk<-as.numeric(unlist(strsplit(as.character(geno[z]), split=" ")))
kk1<-as.numeric(unlist(strsplit(as.character(expr[z]), split=" ")))
d1 <- data.frame(geno1=kk,expr1=kk1)
glm1 <- glm(expr ~ geno, data=d1,family=gaussian(link="log"))
pval <-summary(glm1)$coefficients["geno","Pr(>|t|)"][1]
z<-z+1
}
#results(plot(glm1))
